Magnetic Field Induced Low-Energy Spin Excitations in YBa2Cu408 Measured by 

High Field GS+ ESR 
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We have measured the spin susceptibility of the underdoped high temperature superconductor, 
YBa2Cu4 0g by Gd"^"*" electron spin resonance in single crystals and aligned powders at several 
magnetic fields between 3 and 15.4 T. At low temperatures and high fields, the spin susceptibility of 
the Cu02 planes is enhanced slightly in the -B || c orientation with respect to the B L c orientation. 
The enhancement in an applied field of 15.4 T (~ Q.lbHc2) at 16 K (0.2 Tc) is approximately 10 
percent of the susceptibility measured at Tc- Such a small magnitude suggests that the second 
critical field of superconductivity, Hc2 ~ 100 T, would not suppress the pseudogap. This work 
demonstrates the potential of high field ESR in single crystals for studying high Tc superconductors. 
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The structure of vortex lines in high temperature su- 
perconductors (HTSCs) may shed light on the micro- 
scopic mechanism of superconductivity. The earliest d- 
wave pairing theories implied that the zero temperature 
spin susceptibility should scale with magnetic field as 
(B I H(a)^l'^ in contrast to the linear dependence 

in an s-wave superconductor. A rigorous solution of 
the Bogoliubov-de Gennes equations in the mixed state 
1^ reveals that quasiparticles in d-wave superconductors 
are not bound to vortices, and predict a space aver- 
aged density of states (DOS) iV(0) cx {BjUc^V ^ where 
0.4 < 7 < 1 depending on the temperature and the 
model used . A number of experiments were targeted 
at the field dependence of the DOS in YBa2Cu307_5. 
Scanning tunneling microscopy resolves the DOS at 
and around vortices, but is unsuitable to measure the 
total DOS 0]. Bulk experiments on nearly optimally 
doped YBa2Cu307-5 include heat capacity infrared 
transmission , thermal conductivity , and high field 
muon spin rotation Most of these experiments, to- 
gether with early NMR on planar copper and oxygen [|2j , 
indicated some increase in the DOS with magnetic field 
but were inconclusive about its magnitude. 

Underdoped HTSCs — which also possess a pseudogap 
in the normal state — are of special interest, as it is un- 
clear how the pseudogap is related to the d-wave super- 
conducting gap. An indication whether high magnetic 
fields suppress the pseudogap or not may contribute to 
understanding its nature. Also, most theoretical calcu- 
lations of the quasiparticle spectra in HTSCs depend on 
the assumption that Landau's theory of the Fermi-liquid 
is applicable below Tc, which is still to be justified by ex- 
periments. A recent ^^Cu NMR study on YBa2Cu408 at 
high fields (l3| reports a sizeable field induced DOS for 



B in the (a, 6) plane, implying a large enhancement for 
B II c. In contrast, we show that for B || c the DOS en- 
hancement is small and incompatible with a suppression 
of the pseudogap at i? = H^i ■ 

In this Letter, we report the spin susceptibility of 
YBa2Cu408 obtained by high field electron spin reso- 
nance (ESR) spectroscopy of Gd'^^ in Gd:YBa2Cu408. 
Gd'^+ is a non-perturbing probe of the Cu02 spin polar- 
ization [Q. We search for low-energy spin excitations 
induced by high magnetic fields in the superconducting 
phase. Our goal is to determine the field dependence 
of the anisotropy of the susceptibility, Xs ~ Xs**! from 
the Gd^^^ ESR shift. The anisotropy is a measure of 
the magnetic field induced DOS along c since Hc2 in the 
(a, b) plane is several times larger than in the c direction. 
We measure both oriented powder (OP) and untwinned 
single crystal (SX) samples at several temperatures, fre- 
quencies, and magnetic field orientations. SX data are 
used to determine the zero field splitting (ZFS) parame- 
ters in the spin Hamiltonian of Gd'^+ in Gd:YBa2Cu408 
with precision. This allows us to model powder spec- 
tra and measure shifts in the OP sample in both -B || c 
and i? _L c orientations with high accuracy. Diamagnetic 
effects inhibit SX measurements at low temperatures. 

The experiments were carried out on the high field ESR 
spectrometers in Budapest and Grenoble. In Budapest, 
a quartz oscillator stabilized Gunn-diode was used as a 
mm-wave source at /=75, 150, and 225 GHz. In Greno- 
ble, we used the same setup with a Gunn-oscillator at 95, 
190, and 285 GHz and an optically pumped far infrared 
laser at 349 and 429 GHz, and a sweepable 17 Tesla NMR 
magnet. The frequency of 429 GHz corresponds to 15.4 T 
central Gd'^+ resonance field. The radiation is transmit- 
ted through the sample in an oversized waveguide and 
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no cavity is used. We detect the derivative of the ab- 
sorption with respect to magnetic field. The magnetic 
field is calibrated by a reference sample, BDPA [Q, at 
each sweep. The powder sample Gd:YBa2Cu408 was 
sintered by a standard solid state reaction, with 1% of 
Y substituted by Gd. The powder with a characteris- 
tic grain size of 5 /im was mixed with epoxy resin, and 
aligned in a magnetic field. SX samples of typical dimen- 
sions 1 X 1 X 0.2 mm^^ were grown by self-flux method at 
1100°C in 600 atm of O2. Tc = 80 K in both the SX and 
the OP samples. 

The exchange induced shift of the Gd'^+ ESR lines is 
similar to the ^^Y NMR Knight shift and yields the spin 
susceptibility of the Cu02 planes. We define the Gd'^"'" 
ESR shift, in analogy with the NMR notation as 
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with Bq = hf /gofJ^B, where / is the microwave frequency 
and the arbitrary zero of ^'^K is defined by 50 = 1-9901. 
B!^ is the measured resonance field as deduced from the 
measured positions of the Gd'^"'" lines (after correcting for 
the ZFS, i.e., i?^ would be the resonance field of a similar 
probe with no ZFS). Then 
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K'^iBo, T) - i^o" + ^^Ax?(i?o, T) - B2M, T)/B^, 

(2) 



where Kq — 



(Scd ~ 9o)/9o plays the role of the NMR 
chemical shift. We denote by the "pure" g-factor of 
Gd^+ in Gd:YBa2Cu408 from which the exchange with 
the Cu02 spins is eliminated. (thus Klf) is inde- 

pendent of magnetic field and temperature |l^], and its 
anisotropy is small. The spin susceptibility of the Cu02 
planes is defined by x"(S, T) = M^{B, T)/B where A/f 
is the Cu02 spin magnetization. The Gd'^+ shift due to 
the electronic exchange interaction between Gd'^^ and 
Cu02 is linked to the susceptibility through the con- 
stant '^^A « —15 mole/emu (in analogy with the 
NMR hyperfine constant). We neglect the anisotropy of 
^'^A, since a comparison of the anisotropy of the shifts 
in the normal state YBa2Cu307 and antiferromag- 
netic YBa2Cu306 ||l5| shows that the anisotropy of the 
product '^'^Axs is approximately 5%. B^^^ is the bulk de- 
magnetizing field of the supercurrents in the crystallites. 
In this work we shall estimate the vortex contribution 
to xt high magnetic fields from AB/Bq — {'^'^K'^ — 

c 
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Here '^'^K"' is measured directly (Eq. (|])), and 
is determined from the high field data at high tempera- 
tures. Unlike the *^Y case, B^^^^ is not prohibitively large 
since '^'^A is ten times greater than the ®^Y hyperfine 
constant, ^^A. 

The temperature independent ZFS parameters of the 
nearly pure S state S = 7/2 Gd^+ ion were obtained 
from the SX data. The ZFS of Gd^+ is smaller than the 
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FIG. 1. (a) Gd^+ ESR spectra of a single crystal 
Gd:YBa2Cu408 taken with magnetic field along the three 
crystalline axes at T = 40 K and / = 225 GHz. Neither twin- 
ning nor mosaicity is observed. The field reference BDPA 
is cut out for clarity. Numbers 1, . . . , 7 denote transitions 
1-7/2) ^ 1-5/2), 1+5/2) ^ 1 + 7/2), respectively. Bars 
(I) show simulated line positions. We do not know whether 
ip — 0° belongs to B || a or fe. (b) Positions of the lines cal- 
culated from the spin Hamiltonian of Gd^"*" as a function of 
magnetic field orientation at 225 GHz. 



Zeeman splitting in the fields employed and the spec- 
trum consists of the 7 allowed fine structure transitions, 
whose relative positions depend little on the ESR fre- 
quency but are sensitive to the orientation of the sample 
with respect to magnetic field. Spectra of single crys- 
tals were recorded in the temperature range 30-^70 K. 
Fig. |](a) shows the 225 GHz spectra of an untwinned 
SX with magnetic field along the three crystalline axes 
at 40 K. The principal axes of the spin Hamiltonian of 
the Gd'^+ in its orthorhombic environment coincide with 
the crystalline axes. At 40 K Gd'^^ ESR lines are nar- 
row and about the same width for all transitions showing 
that there are no strains or inhomogeneities in the high 
quality crystals. Although the small microwave penetra- 
tion depth necessitates a few hours of averaging for each 
spectrum, the potential to use Gd^^^ ESR in small sin- 
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gle crystals is clear, e.g., to measure internal fields around 
impurity atoms. The use of resonant cavities renders con- 
ventional X-band ESR spectrometers unsuitable to study 
the superconducting state. Although there are several 
reports on Gd'^^ ESR in perovskites, we do not know 
of any other report on the ESR of HTSC single crys- 
tals below Tc- Measuring the positions of the lines with 
magnetic fields in several orientations and at three fre- 
quencies allows us to fit the ZFS parameters of Gd^+ . 
Details of the spin Hamiltonian in a similar compound, 
YBa2Cu306+£c, were published elsewhere The spin 
Hamiltonian describes the SX spectra well, the differ- 
ence between the calculated and measured line positions 
is always less than the line width. The simulated line 
positions as a function of orientation used for evaluating 
the OP spectra are indicated in Fig. |^(b). 

Once the ZFS parameters are obtained from the SX, 
the Gd3+ shifts, ^'^/^"(B, T), may be deduced with high 
accuracy from the OP spectra. To obtain the increase of 
Xs with magnetic field we use ^d^-a measured at sev- 
eral T and B in both a = c and ab directions as shown in 
Fig. Ij. These data are uncorrected for reversible diamag- 
netic fields, while a small irreversibility in iJ^j^^ at low T 
and low B is eliminated by averaging '^'^K" of spectra 
taken with field swept up and down at the same temper- 
ature At low T and high B the central lines fade 
away, we could therefore measure the shift in 15.4 T reli- 
ably only above 15 K. The dependence of the normalized 
shift on the resonance frequency (Fig. ^) is in a large part 
caused by diamagnetic shielding. However, as shown be- 
low, in the c direction at high fields there is also some 
contribution from the vortex spin susceptibility. 

We assume that only xt ^'"^ ^ measurable field depen- 
dence at the fields employed, and that of xf' may be 
neglected. This is consistent with our main result be- 
low that the field dependence of xl weak, since even if 
there is some field dependence of Xs''; it is significantly 
smaller than that of Xs ) because of the much higher H^^'^ . 
Diamagnetic corrections to B^^ are also smaller than to 
i?^ as Ac ^ Aob holds for the penetration depths. The 
measured field dependence of '^'^K°-^ is weak (Fig. ||(b)), 
and at 15.4 T we neglect the diamagnetic correction. At 
15.4 T the raw ah data are already a good approximation 
of the spin susceptibility. Indeed, the shift at 15.4 T de- 
creases with temperature at low T and the slope is close 
to that expected for a pure d-wave superconductor. In 
the weak coupling limit the low temperature slope of the 
susceptibility of a d-wave superconductor with Tc = 80 K 
corresponds to a.'^'^K shift of lOppm/K p^ ]. 

Diamagnetic corrections are not negligible for lower 
fields and in the c direction, and we estimate the field 
dependence of B"^^^ from classical expressions and recent 
numerical results [p"8|-p2[ . Above the vortex lattice melt- 
ing temperature the magnetic field inhomogeneity due 
to vortices is motionally averaged and S^ia equals the 
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FIG. 2. Normalized GA^+ ESR shifts, ^'^K"[B, T), in the 
OP sample as a function of temperature at several magnetic 
fields, (a) oi = c, (b) a = ab. 



bulk demagnetization, which is roughly proportional to 
Iti{Hc2/B) ||i|-|0|. Therefore one expects a 15-30% re- 
duction in Bj;^ from 8 to 15 T. Our bulk susceptibil- 
ity measurements on the OP sample agree with such a 
reduction. A shift, ABgaddie, in addition to demagneti- 
zation appears when the vortex lattice freezes, because 
the magnetic field distribution is asymmetric and peaked 
at the saddle point, which is different from the aver- 
age field. For zero superconducting coherence length. 



U = 0, AB^ll^^^ « 0.037$o/A2, « 3mT in a triangu- 
lar vortex lattice with Xab — 160 nm |2|]. However, nu- 
merical studies show that already in fields where 
the distance between vortices is much larger than £^ab, 
ASgaddiG is reduced significantly. Corti et al. |Q found 
a 3 mT field distribution due to vortices in a YBa2Cu408 
aligned powder for B || c at 9.4 T below 10 K, and the 
corresponding shift is expected to be a fraction of this 
value. 

Now we estimate the increase of Xs ^ = 15.4 T 
and T — 16 K, the lowest temperature where experi- 
mental uncertainties are small at all fields. Fig. |^ shows 
AB = {^'^K''-^'^K''^)Bn-{K^-K^^)Bo, the anisotropy 
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FIG. 3. Field dependence of the spin susceptibility and 
the diamagnetic correction at 16 K in absolute field scale. The 
solid line is a rough estimate of Bdia{B, 16 K) (see text). The 
error bars represent the maximum error in the measurement. 
Most of the error is caused by the uncertainty in Kq — Kq'', 
therefore it will be strongly correlated between experimental 
points. The brace shows 10% of the Gd'^"'' shift measured in 
15.4 T at Tc. 



approach to explain the microscopic origin of the pseudo- 
gap is that it is associated with superconducting fluctua- 
tions above Tc |2^. In a naive picture, in which the nor- 
mal state spin gap were due to incoherent Cooper pairs, 
breaking of the pairs by a magnetic field would restore 
the spin susceptibility to the value measured at T ^ Tc, 
which is approximately 3 times as much as Xs(?c). Our 
result differs from that of Zheng et al. |l3) since they find 
an enhancement of the spin susceptibility with B in the 
(a, b) plane of similar magnitude as we do for B || c. 

In conclusion we found that at T <C Tc an applied 
field of w 0.15Hc2 enhanced the spin susceptibility of the 
underdoped YBa2Cu408 by only « 3% (« 10%) of its 
normal state susceptibility measured at room tempera- 
ture (at Tc). This suggests that even an apphed field of 
Hc2 would not destroy the pseudogap. 

We are indebted to J. R. Cooper for the static suscepti- 
bility measurements, and I. Tiitto for useful discussions. 
Support from the US-NSF-STCS (DMR-91-20000) and 
the Hungarian state grants OTKA 029150 and AKP 97- 
39-22 is acknowledged. 



of the shifts as a function of field. In AB the anisotropy 
of Kq is taken into account but it is uncorrected for 
reversible diamagnetic effects. The anisotropy of the 
"chemical shift" - K^^ = 500 ± 70 ppm was mea- 
sured from the Bq = 15.4 T data at 0.9 Tc. Above this 
temperature a Korringa-like relaxational broadening of 
the Gd'^"'" lines reduces precision. This value of Kq — K§'' 
is consistent with what is found in the YBa2Cu306+a: 
family ||lj,|l^. Fig. || shows that there is a small field 
dependent contribution to Xs(^o, T) - xfiBo, T) = 
(AB - B^.^{Bo, T))/°'^ABo. AB varies httle with field, 
at higher fields it is constant or increases slightly. In 
Fig. H we illustrate the variation of the reversible diamag- 
netic shift using M^iaiB) = A/dia,o In(i?c2/-B), assuming 
Hc2 = 100 T and that the field dependence of Xs is neg- 
ligible below 3T, i.e., we slightly overestimate Bdia- At 
15.4 T the shift anisotropy is AB = 4.5 ± 1.5 mT, while 
■^dia ~ 2.3 mT. As seen in Fig. ||, the vortex contribu- 
tion to XsiBa = 15.4 T, T = 16K) is of the order of 10% 
of the normal state spin susceptibility at Tc- 

According to the semiclassical scaling relations for 
the spin susceptibility derived by Kopnin and Volovik 
|1|, and the more sophisticated study of Ichioka et al. 
|5|, in a d-wave superconductor Xs oc {B/Hc2)'' with 
7 — 0.4 in the low temperature regime, and 7 w 1 in 
the low field regime. We are close to the crossover [|4j, 
15.4T/77c2 ~ 16K/rc, but whatever 7 is in our case, 
the 10% enhancement implies that the spin susceptibil- 
ity of the underdoped YBa2Cu408 below Tc is restored 
to Xs{Tc) at maximum in a magnetic field of i?c2- Thus 
the pseudogap is little affected by magnetic fields that 
suppress superconductivity. Nowadays the most popular 
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